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Abstract 

The reaction rate of the nuclear fusion d^Li ^Be* {2~^ ,0) is estimated in 
the case where the nuclei are confined to each other by a muon. For the 
description of nuclear transitions, a method which is analogous to the Linear 
Combination of Atomic Orbitals has been used. Using the complex coordinate 
rotation method, we found that a molecular d—^—^Li state exists with energy 
(-20.3084 - i0.0066) eV and |^'m(0)| = 0.44 10"^ fm't. The nuclear wave 
functions needed, were constructed in the form of antisymmetrized products 
of harmonic-oscillator functions for the three-cluster approximation, (dda), 
to the five-body (NNNNa) problem. It was found that the reaction rate A 
is strongly dependent on the energy gap between the d^Li threshold and the 
energy of the final ^Be* resonant state. The value of A obtained by averaging 
over the width of this resonance, is 0.183 lO^'^sec"^. 



PACS numbers: 36.10.-k, 33.90.+h, 36.40.+d 



I. INTRODUCTION 



Low-energy interactions among atomic nuclei are of great interest from pure theoretical 
as well as from practical point of view. From the former, they provide valuable information 
on the nucleon-nucleon interaction and on fine details concerning nuclear structure. By this 
we mean the manifestation of the main symmetries of strong interactions, such as charge 
symmetry, time invariance, etc at very low energies. Another issue here is the importance 
of the small components of nuclear wave functions. Very often these small components play 
a crucial role at super low energies we are interested in. From the practical point of 
view, we mention here the importance of the low-energy nucleus-nucleus data in nuclear 
astrophysics. 



*Permanent address: Joint Institute for Nuclear Research, Dubna, 141980, Russia 



1 



Nucleus- nucleus collision data are not available at energies below the keV region 0. 
Meanwhile nuclear astrophysics deals with nuclear reactions occuring inside the stellar ob- 
jects at these low energies and thus the available scattering data are usually extrapolated 
there. However, in certain systems of light nuclei, we may study nuclear reactions at very 
low energies, namely, reactions between nuclei confined inside muonic molecules. Being in 
a molecular state, nuclei have rather small relative kinetic energy and any possible reaction 
can be treated threshold reaction. 

Nuclei in a molecule are relatively far from each other and the probability of their col- 
lision is very small. However, in some cases this is partly compensated by the large value 
of the nucleus-nucleus reaction cross-section which is due to the existence of a threshold 
resonance. In this paper we consider one such case, namely, the fusion of d and ^Li nuclei, 
which forms the threshold resonant state ^Be*{2~^ ,0). 

Previously, we investigated three such cases, namely the dfi'^Be, pfi^'^Be, and tjJ'He* 
0-0. It was found that the reaction rates in all cases studied were rather high indicating 
that if such molecules are formed, the mixing of states could generate nuclear reactions at 
appreciable rate. 

The paper is organized as follows: In Sec. II we apply an LCAO motivated approach 
to derive formal expressions for the reaction rate; in Sec. HI we briefiy describe how the 
three-body Coulomb problem d — fi — ^Li is treated; Sec. IV contains details concerning the 
employed nuclear wave functions. The potentials used and the procedure of calculating their 
matrix elements are described in Sec. V. Finally, our results are presented and discussed in 
Sec. VI. 

II. REACTION RATE 

Consider the molecular system d — fi — ^Li in its S'-wave ground state. Using the complex 
coordinate rotation method (see forth Sec. HI), we found that a quasistationary state exists 
with energy (—20.3084 — zO. 0066) eV which is in accordance with the early estimates of -19.8 
eV § and -23.8 eV 0. 

This binding energy is negligible on the nuclear energy scale, and beyond the typical 
accuracy with which the nuclear spectra are known. Thus, we can assume that the nuclear 
subsystem, d^Li, of the above molecule, is close to its threshold energy. This threshold 
lies 22.2798 MeV above the ^Be ground state i.e. almost at the same energy where the 
^Be nucleus has an excited state at (22.2 — i0.8)MeV ^ with quantum numbers {J'",T) = 
(2"'",0). This state can be formed also from the corresponding nuclear quantum numbers 
d{l^,0), ^Li{l~^,0), and zero nucleus-nucleus angular momentum. Therefore, the nuclear 
subsystem d %i of the molecule d — — %i can be fused via the reaction 

d + %i — > ^Be + 79.78 keV. (1) 

To calculate the reaction rate of this process we employ a method which is motivated by 
the so-called Linear Combination of Atomic Orbitals (LCAO) approach. We have used 
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this method previously to study the d — jj, — ''Be and p — fi ~ ^^Be 0] systems. 

Following the LCAO approach, we write the total Hamiltonian H , of the nine-body sys- 
tem considered here, in two equivalent forms corresponding to the initial and final clustering 

H = Hi + V\ (2) 

H = H2 + (3) 

where 

H, = Hd + H^L, + Ho + V' + i/^ (4) 

and 

H2 = HsBe. (5) 

are the 'unperturbed' or 'channel' Hamiltonians. In these expressions, H^, Hen, and Hs^e 
are the nuclear hamiltonians of the clusters; Hq is the kinetic energy operator describing the 
nucleus-nucleus relative motion while is part of the strong potentials not included in 
and Hex^i, i.e it represents interactions among nucleons belonging to different clusters; is 
the nucleus-nucleus Coulomb potential and consists of the muon kinetic energy operator 
and its Coulomb interaction with the nuclear particles. 

The 'unperturbed' Hamiltonians Hi and H2 have eigenstates \E'i = A{'^d'^^Li'^moi} and 
^2 = ^8Be*0M , satisfying 

Hi^i = Ei^i, 1 = 1,2 (6) 

and eigenenergies Ei and E2 which are very close to each other. In the above '^moi is the 
three-body wave function of the molecular state; 0^ is the atomic wave function of the 
muon in the Coulomb field of ^Be and '^eii-, and '^s.Be* are the fully antisymmetric wave 
functions of the nuclear subsystems while A is the antisymmetrizer involving nucleons from 
different nuclei. 

In this system, the muon merely acts as a mediator creating an effective force which 
holds the nuclear clusters, d and ^Lz, together at molecular distances. In order to gain a 
qualitative insight into the problem, one can assume that the effective interaction between d 
and ^Li is as schematically shown in Fig. 1. When the clusters are close to each other they 
merge into ^Be* which can be considered as an eigenstate of the deep well. When they are 
separated at large distances, they can occupy the molecular eigenstate of the shallow well. 
The basic feature of this system is that the eigenvalues of these two eigenstates Ei and E2 
are practically the same. As a result, the molecular and the compound nucleus states can 
be mixed. Despite the simplicity of such a two-cluster picture, we expect that the mixing of 
the two states can in reality take place as well. 

The 'perturbations' and mix the 'unperturbed' solutions \I'i and \I'2 and shift 
the eigenenergies Ei and E2. Thus, from the above arguments and following the LCAO 
approach, we look for a solution of the Schrodinger equation 
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(7) 



in the form of a linear combination 



\1/ = Ci^^i + C2^E'2- 



(8) 



Using Eq. (^) in Eq. (0) and projecting with (\[^t|, we obtain the following algebraic linear 
equations for the coefficients Cj [|] 



Ci(i/n-^) + C2(i/i2-^/)=0, 



Ci(i721 - ^/*) + C2(i/22 - ^) = 0, 



(9) 



where Hij = (\E'j|if |\E'j) and I = (\l'i|\E'2), the overlapping integral. The determinant of the 
degenerated system (|^) is equal to zero for the two values of the energy-parameter 

(10) 

[Hn + H22 - {IH21 + I*HuW - 4(1 - |/|2)(ifn^22 - ^^12^^21) }, 

which are the 'perturbed' energy levels corresponding to the 'unperturbed' values Ei and 

E2. 

To find out which of the two levels E^^^ and E^~^ stems from the molecular and which 
from the nuclear state, we modify Eq. (|TD|) by writing 

(11) 



Hu = E, + V,\, 



H12 = Ell + V12, 
if 21 = Ell* + V21, 

'22; 



7^22 = E2 + H2 



where Vfj = and H^2 = (^2|i^''|^2). Noting that (1 - |/ 

g = |/p + l/l"^ + |/|^ + ... is a small quantity, we get 

E(-) = E2 + A2, 

where 



2\-l 



(12) 
(13) 
(14) 

I + g, where 



(15) 



Ai = Vi\ + gHii + (1 + (?) (VT+s - l)iEi -E2 + Vi\ - ^2) ~ H^H: 



21, 



and 



A2 = H^2 + ^?^22 ~{i + g) {VT+s - l)iEi -E2 + Vi\ - H^2) + MIH21] 



\I\'HiiH22 + \H2i\' - {Hii + H22)Ft{IH2i) + {Im{IH2i)f 



{El -E2 + Vi\ - ^2^ 



It is clear that is the 'perturbed' molecular energy level, and Ai is the shift caused by 
the strong interaction between the nuclei, while A2 is a negligible shift of the nuclear level 
caused by the muon cloud surrounding the ^Be* nucleus. 

From the above expressions for Ai and A2 it is seen that our approach differs from the 
simple perturbation theory where Ai ^ V^i and A2 ~ The LCAO wave function is, 
therefore, written 

^(+) = Ci+^^i + C^+^^s, (16) 
where C['^^ and c!^'^^ are solutions of (P) corresponding to the energy E = E^^K 
The probability P for the transition ([|) is obtained as usual by the projection 



P 



(^2!^^+^) . (17) 



The P can be interpreted as the transition probability through the potential barrier. It 
is noted that the total wave function '^^'^^ exp{—iE^^H/h) oscillates at the barrier with 
frequency z/ = \E^'^^\/2iTh. Hence, in order to obtain the reaction rate A, we must multiply 
the probability P by u. Solving the algebraic system (^, we finally get 

2nh |/|2 + //* + /*/ + !' ^ ^ 

where 

Therefore to calculate the reaction rate we need the matrix elements V(y-, V^i2) -^^22; ^'^'^ 
the overlapping integral / = ^^^\^2}. This in turn requires the construction of the 'unper- 
turbed' molecular [d — fi — ^Li) and nuclear i^Be*) wave functions and ^2 respectively. 



III. MOLECULAR STATE 

In the initial state, the nuclei d and ^Li are very far apart and held by the muon at 
molecular distances (in the shallow well of Fig. 1). Since this separation is much greater 
than the nuclear sizes, we consider them as point-like charged particles interacting via a pure 
Coulomb potential. In order to solve the resulting three-body {dix%i) Coulomb problem, we 
employ the complex coordinate rotation method which has been successfully used in atomic 



physics for calculating positions and widths of resonant states ( see, for example, Ref. |T0 
and references therein ). 

The main idea of this method consists in the complex rotation of the spatial coordinates 

r ^ fexp(^), e = er + iOi, (20) 
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which makes the resonance state wave function square integrable. This, in turn, permits 
the use of variational methods for their location since the rotation leaves the resonance pole 
positions intact. Another advantage of this approach is that the search for resonances as 



eigenvalues in the complex plane has a rigorous mathematical foundation [p!3| ,p!4 



Muonic molecular systems have extremely narrow resonances, with widths which are 
difficult to calculate, because the accuracy of the numerical calculation for the complex 
energy has to be at least higher than the width magnitude. A variational method based 
on a random-tempered Slater-type exponential expansion exhibited high convergency and 
accuracy in calculations of muonic molecular bound states Here we demonstrate 

that this method is also successful in the case of dfi^Li system. 

The system of interest consists of three particles, a negative muon of mass and two 
nuclei of masses Ma and M),, where a, hereafter, stands for a lithium nucleus and b for the 
deuteron. The Hamiltonian (in muonic atomic units e = h = = 1), after separating the 
center of mass motion, can be written as 

where fa and fb are the vectors towards the muon from the two nuclei, R denotes the distance 
between the nuclei, and = m^Mi/{m^ + Mj) are the reduced masses of the respective 
muonic atoms (z = a,b). This set of coordinates, fa, fb, and R, is more appropriate to handle 
the molecular three-body problem than the Jacobi coordinates {R, f^) commonly used in 
few-body problems. 

In the complex coordinate rotation theory the resonant state is defined as the solution 
of the eigenvalue problem 

{H^om - E)m^oi = 0, Hmom = U{e)HmoiU-\e), (22) 

for the Hamiltonian Hmoi{0) analytically continued into the complex plane of dilation pa- 
rameter 9 of the transformation (pO|). The corresponding transformation of the molecular 
states in the Hilbert space is defined via 

[U{e)F] (f) = e=^'/'F(eV). (23) 

Such a transformation has a great computational advantage for systems with Coulomb 
interactions. The kinetic and potential parts scale as exp {—29) and exp (—6*), respectively, 
and the Hamiltonian can be written as 

H^oi = Te-^' + Ve-'. (24) 

The variational wave function for the state with total orbital angular momentum L = 
is expanded in the form 

N 

"^moiiR, fix) = ^Ci exp(-aira - I3in - liR), (25) 

i=l 
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where the vector directed from the center of mass of the nuclear subsystem to the muon, 
is the hnear combination 

Mafa + Mbfb 



r 



Ma + Mb 

of fa and fi,. The nonhnear parameters are chosen by a pseudorandom algorithm, previously 
used by Thakkar and Smith |jl5|, and are generated by the formulae 



ai= A{V2[i{i + l)/2]), 

A=B(v^[z(z + l)/2]), (26) 
7.= C(v^[^(^ + l)/2]), 

where (. . .) denotes the fractional part of a number. The tempering parameters A, B, and 
C can be taken more or less arbitrarily. A resonable requirement is that these parameters 
must be chosen to reflect the geometry of the molecule. 

Once the nonlinear parameters Oj, and 7^ are fixed, the coefficients q can be found 
by solving the complex nonhermitian matrix eigenvalue problem 

Acfc = E,,Bcfc, (27) 

where 

A = {'^mol\HmoliO)\^mol), B = (^'^oZ | ^^oz) • 

The corresponding discrete complex eigenvalues have the form 

Ek = Er- tTm/2 (28) 
where Er gives the position and the width of the resonance. 

It should be emphasized that the exact eigenvalues do not depend on the rotating param- 
eter 6. However, employing the ansatz (^5]), we deal with approximate eigenvalues which, 
therefore, have such dependence. On the other hand, in the area near the spectral points 
(including resonances) the dependence of eigenvalues on any parameter should be weak since 
the corresponding functional must be stable. Since the variational parameters of the expan- 
sion (p5D are fixed by (^6D and (p7|) , then the only parameters which remain free are the real 
and imaginary parts of dilation angle 6. Therefore, the resonance positions and widths are 
deduced from the condition that a discrete complex eigenvalue is stabilized, i.e 

dEk , , 

with respect to variations of the complex dilation parameter Q. 

To solve Eq. (^7[) numerically, the inverse iteration method adapted to symmetric com- 
plex matrices has been employed. This method is very effective, as far as computational 
effort is concerned, and stable to round-off errors. This enables us to extend the number of 
the expansion functions up to 1400 terms. 
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In our approach, the full three-body wave function '^moi is not needed since in actual 
calculations of the fusion rate, only the effective two-body wave function describing 
the nucleus-nucleus relative motion is required. Moreover, will appear only in the integral 
expressions involving either the short-range strong potentials or the nuclear wave-functions 
localized within a volume extending to few fm. Hence, such integrals take into account only 
a small area around the R = point. The boundary value "^mi^) is obtained from "^moi 
by noting that after the complex rotation, the resonance wave function is square integrable 
and thus, the density p^^\^) of ^^e probability for the nuclei to be at vanishing internuclear 
separation, i? = 0, is given by 

mol I ^mol) 

where "^Inoi is the solution of the adjoint equation 

^l,,H^a{e)=^LiE*. (31) 
Finally, the rotated value of p^^\^) can be easily transformed to 6 = via 

pff(0)=e-^Vff(0). (32) 
Therefore, the boundary value \E'm(0) is 



1*^0)1 = vW(0)- 



Our calculations were performed in a 32 digits arithmetics. The values of Me^^j , Md and 
are (in units of electron mass) 10961.9, 3670.481 and 206.7686, respectively, while the values 
of the tempering parameters used are A=1.8, B=1.8, and C=2.8. 

In Fig. 2 we present the results obtained for an eigenvalue as a function of the complex 
dilation parameter and with variational wave function (^3]) having a total number of terms 

= 900. The curves correspond to different values of the parameter a = e^''. Each 
curve connects points on the rotational path for the fixed value of a and different values 
of the rotational angle 6i = 0,0.02, . . . ,0.2 in steps of A6i = 0.02. It is clearly seen that 
the rotational paths form the stabilized eigenvalue of (—20.3084 — z0.0066)eV at 6i ~ 0.1, 
the accuracy being of the order of 10~^eV. The F^ corresponds to the predissociation rate 
of 1.810^%ec~^. The boundary value of the nucleus-nucleus wave-function obtained by the 
above procedure is 

1^^(0)1 = 0.44 10"^ fm-^/^ 



IV. NUCLEAR MODELS 

In order to reduce the eight-body nuclear system, d^Li, in the initial channel to man- 
ageable proportions, ^Li is considered to be a ((ia)-bound state system. In this way we have 
in the initial state a three cluster system, namely, the dda one. It is known that such a 
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two-body assumption for the %i ground state, is sufficiently accurate |]T8[. This model is also 
consistent with the treatment of the final state of the transition d^Li ^Be*, since when 
^Be nucleus is excited at ~ 20MeV, then only one of the two a-clusters in it is destroyed 
IITSl forming two deuterons. This is schematically depicted in Fig. 3 together with two sets 
of possible Jacobi coordinates. 

The complexity of the problem, especially the required multidimensional integrals (up 
to dimension of 10), necessitated the use of bound state wave functions with simple radial 
dependence which, however, are widely employed in the literature. These wave functions 
describe reasonably well the nuclear sizes which are crucial in our model. We shall discuss 
them below in more details. 



A. The deuteron wave function 

This is taken to be of Gaussian form 

Nd ( d 



v|/,(f) = -/=exp ( --r^ 1 xV, (33) 



with 



and 



where and t]"' are the spin symmetric and isospin antisymmetric functions. The mean 
square radius is fixed by (r^) = (1.956 fm)^ [O]. 



B. The ^Li wave function 

This is constructed from the product of and the harmonic oscillator 2s- wave function, 
describing the motion of the a particle with respect to the center of mass of the deuteron. 
It can be shown that, in the six-nucleon oscillator model of ^Li-ground state, the degrees of 
freedom associated with the variable r2 (see Fig. 3), has two quanta of excitation [IE\. This 



is a consequence of the Pauli principle. Therefore, by using the 2s-function we, indirectly, 
take into account effects of the Pauli principle despite the fact that the a-particle is treated 
as elementary. Therefore the following ansatz was used 

^^Liin, fa) = *d(fl)02s(f2), (34) 

with the 2s-function having the usual form 



<P2s{r2) = -j^^ 4 - - exp --^ . (35) 
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The oscillator parameter a, 




is chosen to reproduce the experimental mean square radius of ^Li, (re^.) = (2.54 fm)^ [|^. 



C. The ®Se*(2+,0) wave function 

As mentioned earlier, the wave function of ^Be*{2~^ ,0) state, is assumed to be that 
of the five-body system {nnppa). It is well known that the ground state of ^Be can be 
described reasonably well by the two a-cluster model [|l8l. Using the eight-body harmonic 
oscillator model, however, this ground state is a (ls)^(lp)^-configuration having four quanta 
of excitation. On the other hand, in the a — a model all four quanta must be attributed to 
the a — a relative motion because each of the two a-particles is in the (Is) ^-configuration. 
Hence, for the ground state of ®Be the function Al"^ ai"^ a2(p3s} is expected to be a good 
approximation. Here 



1 2V30 



A i i 

47r 7r4C2 



1 1 /r\2 1 /r^ 4 
4 ~ 3 Vcj ^15 Vc 



^2 



is the normalized harmonic-oscillator function having two nodes and describing the relative 
a — a motion with four quanta of excitation. The oscillator parameter c for this function, 

c^ = ^((rL)-(0), 

is chosen to reproduce the experimental size (r^g) = (2.39 fm)^ of Berillium nucleus pO 
For (r^) we have used the value (1.671 fm)^ . 



At 22.2 MeV above the ground state, the whole excitation energy is used to break one 



of the a-clusters [0. It is therefore natural to assume that the other a particle remains 
in the same 3s-state with respect to the center of mass of the destroyed cluster and bears 
the whole excitation. Consequently the excited a cluster should have quantum numbers 
corresponding to the (2+, 0) state. The positive parity is obtained if we have two quanta of 
excitation in the oscillator model. The simplest configuration in this case is the dd-system 
with relative motion wave function (^2s- The oscillator parameter in this case is adjusted in 
order to get the correct value of the excitation energy. 

Therefore, the wave function of the resonant state '^_Se*(2"'', 0) is 

^8Be*(n,r4,f5,r6) = Ar203^(f6)^{^d(fi)*rf(f4)02s(r5)^(12)(34)}, (36) 

where 

e(i2)(34) = r^?2%"4 E(l^l^'|2, a + a')xl2i^)xU^') (37) 
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is the four-nucleon spin-isospin wave function constructed from the antisymmetric isospin 
functions rjf^ and the symmetric spin functions xtj of the ij-nucleon pair and A/2 is the 
normahzation constant; A is the antisymmetrizer consisting of the following permutations 

A = (12)(34) - (32)(14) - (42)(31) - (13)(24) - (14)(32) + (34)(12). (38) 

It is known that the antisymmetrization significantly reduces the differences among the var- 
ious cluster representations of nuclei [|1^. Moreover, in the harmonic oscillator model all 



possible cluster representations, after being antisymmetrized, become equivalent. This give 
us reasons to believe that our wave function for ^Be*, constructed in the cluster representa- 
tion dda, is a reliable one. 



D. The final and initial states 



The wave function "^SBe* appearing in the final state \E'2 can be obtained from (PB]). After 
carrying out the normalization and antisymmetrization it is written as follows 



?e* iri,r4,r5,r6j 



144 3,3(1 - 2db' 



3 ^d2b^- — — 

V2 (1 + 2rf62)5_ 

hsire) \ ^d(ri)^'d(f4)02s(r5)^(l2)(34) 



-*<i(y - ^5 + + ^5 + y)02s( )^(32)(14) 

-^d{ — r5 - -Jj'^di— + - — )02s( )^(13)(24)|- 

The initial state wave function \E'i = A{^d^<^Li'^moi} can be similarly antisymmetrized with 
the help of the operator given in Eq. (B3). Thus, after some algebra, we obtain 



^i(fi,f2,r3,r4,f7,) = 



T2 5r2 r3 

^2s{r2)'^mol{r3,rf,) + (p2s{— - n!)^moi( TT, ^m) 



6 



^(12)(34) 



^ ,fi 2r2 ^ rA^ ^ ,ri 2r2 ^ Ta, 
dv 2 3 d 2 ^ 2 3 2 ' 



w ^1 ^ 2r2 
52^( — r + 



2 + ^)^mol[ 3 9 + 3 + 3'^/^J 



4 3 

fi 2f2 rs u fi 4f2 fs f4 ^ 

dl 2 3 2' 2 3 2' 



^(14) (32) 



!>2sl 



fi 2r^ _ _ _ 1^2 _ ^ ^ x 

•4+3 2 4^^™"'^ 3 9+3 3'^/^^ 



fi 2f2 fg f4 fi 4f2 fg r4 ^ 

+ 02.( J + ^ - y + j)^^oKy - ^ + y + Y'-.) 

The normalization constant A/i is obtained numerically. 



^(13)(24) \ ■ 
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V. POTENTIALS AND MATRIX ELEMENTS 



For the description of the strong interactions among the particles of the five-body 
[NNNNa) system we need the nucleon-nucleon, v^^ , and nucleon-a, f^^, potentials 
(m,n = 1,2,3,4). 

Since and \l/2 are constructed from S'-wave components only, the total angular mo- 
mentum J = 2 of the (NNNNa) system can be formed from the spin-momenta of the 
four nucleons. This is possible when the spins are aligned, i.e. when the total spin of any 
A^A^-pair equals 1. Thus, as a nucleon-nucleon potential, we employ the Malfliet-Tjon triplet 
potential [p^ , 

v^^{r) = V,^-^-V,^-^, (39) 
with Vi = 1438.72 MeVfm, V2 = 626.885 MeVfm, ai = 3.11 fm-\ and as = 1.55 fm'^ 



For the Na potential, we choose the one proposed in Ref. [p^ ] 

v^^^r) = -Voexp{~wr^) (40) 
where Vq = 55.774 MeV and w = 0.292 fm'^. 

For the matrix elements we have 

^ mv^/ + + + <^ + vl" + (41) 

Due to the symmetry properties of \E'i and with respect to nucleon permutations, we 
may write 

V^; = 4(vl/,K/|v[/^.) + 2{^,\vl-\^,). (42) 

After antisymmetrization, the wave functions \E'i and \E'2 acquire components with nonzero 
orbital angular momenta. However, in the integrals (|42|) the S'-wave potential operators 
which implicitly include the projection operators 7^(^14 = 0) and V{i4,a = 0) which retain 
only S'-waves along the vectors 



X = h H , 

2 3 ^ 2 ' 



r2 ^ 
l/=-y + r3 + -, 



(43) 



towards the first nucleon and the a-particle from the position of the fourth nucleon respec- 
tively (see Fig. 3 ). In order to perform this S'-wave projection we make the additional 
integration 



V{iu = 0)^. = ^ / = "^"^ 
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Then the matrix elements are calculated as follows 



oo 

= (477)2 J dr^d^^r^xf J dr2df4dr^^lv^/'^,^, (44) 



oo 

= (47r)2 J dndyinyf J dr,dndf,%vl"<^>,y. (45) 



For the overlapping integral I we have 

oo 

I = An J drir^ J dr2dridrdr^^!\{fi,r2, ns, r^,)'^SBe* (n, ^4, + ^3 



2r2 , ^ 2r2 rs, 
3~ ~ 2"' 



(46) 



where ri is taken along the z-axis of the reference frame for the vectors r2, ra, and r4. The 
dot products among these vectors, appearing in the integrand have the general structure 

rm-rn = rmTn [siu Om siu On COs(v9m - V?n) + COS 6^ COS 6^] , 

where m, n = 1, 2, 3, 4 and ^m, are the spherical angles of fvi in that frame. The combi- 
nation '^\'^SBe* contains dot products among the spin-isospin functions 

6 = ^(12)(34), 6 = ^(14)(32), and ^3 = ^{13)(24)- 

It is easily seen that {ik\ii) = 1/2 for all k and /. 



The above multidimensional integrals (0), (45), and (|46|) involve an integration over 



the muon variable (see Fig. 3). In order to reduce their dimensions, we exploit the fact 
that nuclear wave functions and potentials are localized within a small volume of the size 
which is not lager than few fm, while the molecular wave function '$rnoi{R,'^fi) is practically 
constant when R is within a region of ~ 100 fm from the point R = 0. 

Thus, in these integrals we can replace 'ifrnoi{R,r^) by ^l'moz(0, r^). On the other hand, 
^mo^O' ^fj.) is the so-called united atom limit of the three-body molecular wave function and 
describes the motion of the muon in the joint Coulomb field of the two nuclei with zero 
separation between them. Hence, in this limit we have 

^moi{R,r,) -^Jmm,{r,) , (47) 

rC — >U 

where (pfiivfj.) is the atomic wave function, and \l/m(0) is the limit [R —>■ 0) of the wave 
function describing the nucleus-nucleus relative motion inside the molecule ( the calculated 
value of ^'^(O) was given in Sec. Ill ). 

Therefore, by substituting the factorized expression (^) into the above integrals, the 
integration over the muon variable is eliminated. The most formidable of the remaining 
integrals is one for the overlapping I which still has 10 dimensions. 
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Finally, we have to find the matrix element of the Hamiltonian ^^3, ^^4) de- 

scribing the muon motion in the Coulomb field generated by the nucleons. This Hamiltonian 
depends on the nucleon variables parametrically. Thus we have 

and again, exploiting the smallness of the nuclear size, we use the reliable approximation 

= (0^1/7^(0,0,0,0)10^). 

The last matrix element is just the eigenenergy of the muonic atom corresponding to the 
eigenstate |0^) with the quantum numbers {nim). This eigenenergy is the one of a hydrogen- 
like atom and is expressed via the grand orbital quantum number n as —{A/nYRy where 
Ry = m^e^/2h'^ = 2813.25 eV is the Rydberg constant for muon. 

A question then arizes: which of the atomic states should be used here? The answer 
can be found with the help of the Born-Oppenheimer model. Indeed, in this model for the 
hydrogen- lithium system |T^, only the so-called Sda-term provides the shallow well shown 
in Fig. 1. The solutions of the two-center Coulomb problem are usually classified by the 
quantum numbers [nim) of the united atom, {R 0), which is formed when the muon 
moves in the Coulomb field of the total charge of the nuclei. Thus, the symbol 3da means 
that, when the nuclei approach each other, the muon is in the atomic state with grand 
orbital quantum number n = 3, angular momentum i = 2, and m = 0. Moreover, when 
R ^ 00 the 3c/cr-term correspond to free motion of lithium nucleus and (i/i-atom in its 
ground state. Hence, as the zero point of the energy scale we should choose the ground state 
energy, {—Ry), of (i/i-atom. Therefore 

m2^-[l)' Ry + Ry = -lRy- 

We epmhasize that the Born-Oppenheimer model was used here only to find the quantum 
numbers of the atoms which are formed in the limits R ^ and R ^ 00. 



VI. RESULTS AND DISCUSSION 

The numerical evaluation of the multidimensional integrals for the potential matrix el- 
ements and the overlapping integral, were performed with the help of the Haar-function 
expansion method The results thus obtained are 



= -0.14524 10-2 eV, 
= = -6.2484 eV, 
H^^ = -5001.3 eV, 
I = 0.42501 10"^ 

The use of these matrix elements in the formulae of Sec. II along with 

El = 22.2798 MeV and E2 = 22.2 MeV 
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gives for the molecular level displacement Ai, transition probability P, and reaction rate A 
the following results: 



Ai = -0.963 10"^ eV, 

P = 0.606 10-^ (48) 
A = 0.298 10® sec-^ 

At this point it is important to note that the compound nuclear state ^Be*{2~^, 0) is a resonant 
one with energy, (22.2 — iO.8) MeV, having a rather large width F = 1.6 MeV. This width 
is much greater than the assumed difference, 

5 = Ei-E2 = 79800 eV, 

between the molecular and compound nucleus levels. We further note that a resonance can 



be excited within a wide energy interval around the central resonant energy [^, and the 
relative probability of its excitation is defined by the Breit-Wigner factor 



(E-E2)2 + rV4' 



Therefore, in order to obtain the correct values of Ai,P, and A, one should calculate them 
with different values of 6 = Ei — E2 and take the average value with W{E) as the weight 
factor. 



The dependence of Ai, P, and A on the energy level difference 6 is shown in Table |I[ It 
is seen that the reaction rate A significantly increases when 6 —>■ 0, and attains a very high 
value 0.246 10^^ sec~^ when the levels Ei and E2 coincide. To include these contributions 
we note that the nuclear transition (|l|) can only take place if Ei > E2. Further, since 5 
cannot be greater than Ei, the interval of 6 for the averaging is [0,i?i]. Thus, the average 
value of the reaction rate is obtained via 

fdSX{S)WiEi~5) 
(A) = . (49) 

/ dSW{Ei - 6) 


The averages (Ai) and (P) , can be similarly defined. The results which take into account 
the spreading of the final resonant state are 

(Ai) = -0.116 10-2 eV, 
(P) = 0.373 10"^ 
(A) = 0.183 10^° sec-^ 

Our conclusions can be summarized as follows. We employed the LCAO-motivated ap- 
proach to calculate the transition probability and reaction rate for the nuclear fusion inside 
the muonic d — fi — ^Li molecule by using the fact that ^Be has an excited state near the 
d^Li threshold energy. The wave functions of the nuclear subsystems in these calculations 
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were constructed from antisymmetrized products of the harmonic oscillator functions in a 
three cluster (dda) approximation to the five-body model (NNNNa) while the three-body 
molecular problem was treated in the framework of the complex coordinate rotation method. 

Within this model, the reaction rate strongly depends on the gap between the initial and 
final state energies. Since, however, the final state in our case is a resonance, we consider 
it reasonable to average the results over the permissible (by the reaction) energies covered 
by the resonance. The average results turned out to be two orders of magnitude higher as 
compared to those obtained using the center of the resonance. This makes the reaction (1) 
quite attractive for further theoretical as well as experimental investigations. 



Acknowledgements 

Financial support from the University of South Africa and the Joint Institute for Nuclear Re- 
search, Dubna, is greately acknowledged. This work was also partly supported by the Grants 
RFB300 and NSCTOOO of the International Science Foundation and Russian Government, 
and by the NATO Collaborative Research Grant N930102. 



16 



REFERENCES 



[1] C. Rolfs, in Nuclear Astrophysics, Research Reports in Physics, p. 162, Springer Verlag 
(1988). 

[2] V. B. Belyaev at al. Communication of JINR, D15-92-324, Dubna, 1992. 
[3] V. B. Belyaev, H. Fiedeldey, and S. A. Sofianos, Sov. J. Nucl. Phys., 56 (7), 66 (1993). 
[4] V. B. Belyaev, S. A. Rakityanski, H. Fiedeldey, and S. A. Sofianos, Phys. Rev., A50 
(1), 305 (1994). 

[5] V. B. Belyaev, V. I. Korobov, and S. A. Rakityansky, Few-Body Systems, 17, 243 (1994) 
[6] A. V. Kravtsov, N. P. Popov, and G. E. Solyakin, Sov. J. Nucl. Phys., 35 (6), 876(1982). 
[7] V. B. Belyaev et al. • • •, Few-Body Systems, SuppL, 6, 332 (1992). 

[8] F. Ajzenberg-Selove, Nucl. Phys., A490, 1(1988). 

[9] J. C. Slater , Electronic Structure of Molecules, McGraw-HU Book Company, 1963. 
[10] Y.K. Ho, Phys. Reports, 99, 1(1983). 

[11] A. M. Frolov, V. D. Efros, Pis'ma Zh. Eksp. Teor. Fiz. 39, 449(1984) (JETP Lett. 39 
(1984) 545); A. M. Frolov, Yad. Fiz., 44, 1367(1986) (Sov. J. Nucl. Phys., 44 888(1986)). 
[12] S.A Alexander and H.J. Monkhorst, Phys. Rev. A, 38, 26(1988). 

[13] J. Aguilar and J. M. Combes, Commun. Math. Phys. 22, 269(1971); E. Balslev and J. 

M. Combes, ihid. 22, 280(1971). 
[14] B. Simon, Commun. Math. Phys. 27, 1(1972). 
[15] A. J. Thakkar and V. H. Smith, Jr., Phys. Rev. A 15, 1(1977). 
[16] E. Brandas and P. Froelich, Phys. Rev. A, 16, 2207(1977). 

[17] I. V. Komarov, L. I. Ponomarev, and S. Yu. Slavyanov, Spheroidal and Coulomb 

Spheroidal Functions, 'Nauka', Moscow, 1976. 
[18] K. Wildermuth and Y. C. Tang, A Unified Theory of the Nucleus, Vieweg, Braunschweig, 

1977. 

[19] S. Klarsfeld, J. Martorell, and D. W. L. Sprung, J. Phys. G: Nucl. Phys., 10, 165(1984). 

[20] I. Tanihata et al., Phys. Rev. Lett., 55, 2676(1985). 

[21] D. R. Tilley, H. R. Weller, and G.M. Hale, Nucl Phys. A 541, 1(1992). 

[22] B. A. Malfliet and J. A. Tjon , Ann. Phys., 61, 425 (1970). 

[23] S. B. Dubovchenko and A. V. Dzhazairov-Kakhramanov , Sov. J. Nucl. Phys., 51, 

971(1990). 

[24] I. M. Sobol, Multidimensional Integration Formulae and the Haar Functions, Nauka, 
Moscow, 1969. 

[25] V. L Kukulin, V. M. Krasnopolsky,and J. Horacek, Theory of Resonances, Kluwer Aca- 
demic Pubhshers, Dortrecht, 1988. 



17 



TABLES 



(eV) 


Ai (eV) 




A (sec) 





-6.25 


0.499 


0.246 10^° 


1 


5.77 


0.460 


0.226 10^° 


10 


3.00 


0.188 


0.921 10^^ 


10^ 


0.387 


0.386 10-2 


0.190 10^^ 


103 


0.376 10-1 


0.390 10-"^ 


0.192 10^2 


10^ 


0.245 10-2 


0.390 10-6 


0.191 10^° 


79800 


-0.963 10-3 


0.606 10-» 


0.298 10^ 


10^ 


-0.106 10-2 


0.385 10-^ 


0.189 10^ 


10^ 


-0.141 10-2 


0.339 10-1° 


0.167 10^ 



TABLE I. The dependence of the molecnlar level displacement Ai, transition probability P, 
and reaction rate A on the energy difference 6 between the molecular Ei and compound-nucleus 
levels. 
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FIGURES 




FIG. 1. Typical structure of the effective two-body interaction potential acting between the 
nuclei confined within a molecule. 



19 



ImE (IQ-W) 



0.00 



-0.40 - 



-0.80 - 
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FIG. 2. Complex eigenvalue trajectories corresponding to variations of the rotation angle 9i 
and for different values of dilatation parameter a = e^*" . 
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FIG. 3. Two alternative sets of Jacobi vectors {n, r2, ra, r4, f^} and {n, f4, rs, re, r^ii} describ- 
ing space configuration of the six-body system {NNNNa/j,). The x and y are auxihary vectors 
defined by Eq. (43). 
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